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We reexamine basic aspects of a nonequilibrium steady state in the Kondo problem for a quantum 
dot under a bias voltage (eV = — fiR) using a reduced density matrix (RDM). It is obtained 
in the Fock space by integrating out one of the two conduction channels, which can be separated 
from the dot. The remaining subspace is described by a single-channel Anderson model, and the 
statistical distribution is determined by the RDM. At zero temperature, the noninteracting RDM 
is constructed as the mixed states that show a close similarity to the high-temperature distribution 
in equilibrium. Specifically, when the system has an inversion symmetry, the one-particle states in 
an energy region between the two chemical potentials (/i_r < e < are occupied, or unoccupied, 
completely at random with an equal weight. The Coulomb interaction preserves these aspects, and 
the correlation functions can be expressed in a Lehmann-representation form with the mixed-state 
distribution. Furthermore, the statistical weight shows a broad distribution in the many-particle 
Hilbert space, spreading out in a wide energy region over the bias energy eV. 

PACS numbers: 73,63.-b, 72.10.Bg, 73. 21. La 



I. INTRODUCTION 

The Kondo effect in quantum dots has been an ac- 
tive research field over a decade So far, the equi- 
librium and linear-response properties of a single quan- 
tum dot have been understood well, basically, based 
on the knowledge of the Kondo physics in dilute mag- 
netic alloys 2- The nonlinear transport under a finite bias 
voltage V however, is still not fully understood, al- 
though a number of previous work a 10 ! 11 ' 12 ! 13 ' 14 ! 15 and 
recent one o 16 i 17 ' 18 have contributed to the continuous 
progress. One of the complications in the nonequilib- 
rium system is that the density matrix is not a simple 
function of the Hamiltonian H . In thermal equilibrium, 
the density matrix has a universal form p eq oc e~@ H . On 
the other hand, the nonequilibrium density matrix is not 
unique, even for the steady states. 

The Keldysh formalism has widely been used for con- 
structing the nonequilibrium density matrixi 19 : 20 It has 
traditionally been described with a perturbative Green's 
function approach. However, in order to study compre- 
hensively the nonequilibrium properties of strongly cor- 
related electron systems, accurate nonperturbative ap- 
proaches are also required. For instance, in the ther- 
mal equilibrium the approaches such as numerical renor- 
malization group (NRG)2i and quantum Monte Carfo22 
methods played important roles to clarify the Kondo 
physics on the whole energy scale. On the other hand, we 
still have not had systematic ways to handle the nonequi- 
librium density matrix. It makes the construction of non- 
perturbative approaches difficult. Therefore, further con- 
siderations about the basic formulation itself seem to be 
still needed. 

In this paper, we describe a reformulation of a nonequi- 
librium steady state, using a reduced density matrix 
(RDM). It is obtained in the Fock space by integrating 
out one of the two channel degrees of freedom, which 
can be separated out from the dot by a reconstruction 



of the linear combination of the channels. The rest of 
the Hilbert space consists of the dot and the remaining 
channel, which are coupled with each other by a tunnel- 
ing matrix element. The RDM takes a simple form in 
the noninteracting case, as shown in Eq. (|33[) . It enables 
us to treat the nonequilibrium averages using a single- 
channel Anderson model. The bias voltage appears in 
the subspace through an energy-dependent parameter T e 
given in Eq. , which has a close relation to the usual 
nonequilibrium distribution function f c s(e). Specifically 
for the symmetric coupling Tl = ^r, the electrons dis- 
tribute at zero temperature to the one-particle states be- 
tween the two chemical potentials /i£ and completely 
at random. The Coulomb interaction is taken into ac- 
count adiabatically, and then the nonequilibrium aver- 
ages can be expressed in form of Eq. (|48|) . in terms of 
the interacting eigenstates and the RDM. This expres- 
sion could be used for Hamiltonian-based nonperturba- 
tive calculations. It is also deduced that the nonequilib- 
rium statistical weight spreads out in a rather wide en- 
ergy region of the many-particle Hilbert space over the 
bias energy eV. 

The model is described in Sec. HU The derivation of 
the RDM is given in Sec. IIIIl The mixed-state aspects 
of the RDM are examined for T L = T R in Sec. [TV] The 
nonequilibrium distribution in the many-particle Hilbert 
space is discussed in Sec. [V] Summary is given in Sec. 
IVII In the Appendixes, details of some calculations are 
provided. 



II. FORMULATION 

A. Scattering states 

We start with a single Anderson impurity connected to 
two leads on the left (L) and right (R) . The Hamiltonian 
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is given by H = Hq + Hjj with Hq = H d + H c + Ht, 
Hd=[e d + -)n d , Hu = -[n d -lj , (1) 

H c = Yl E d££C t.ua c e,ua , ( 2 ) 
u=L,R <J J 

h t = x; $>„( 4 a, • w 

u=L,R a 

Here, d\ creates an electron with spin a in the dot, 
ii d = J2 a d \r d cn £ d is the onsite potential, and v v is a tun- 
neling matrix element between the dot and lead v. The 
conduction electrons are normalized as {c £ va , c\, — 
Suv' $aa' 5(e — e'), and C vo . = J de^J pie) c e va with a con- 
dition J de p{e) — 1; the same one-particle density of 
states p{e) is assumed for both of the leads. We are using 
explicitly the continuous form of the conduction bands 21 
in order to make the scattering states well defined. The 
Fermi level at equilibrium is taken to be p = 0. The two 
different chemical potentials, p^ = ~Pr = eV/2, are in- 
troduced usually into the isolated leads for Ht = 0<2£ Al- 
ternatively, as described by Hershfield^ the same steady 
state can be constructed from Hq including Ht in the 
initial condition, and it has been applied to some special 
models ! 12 i 24 We use the later formulation in the following 
to set up the nonequilibrium density matrix, and assign 
Pl and pu to the scattering states incident from left and 
right, respectively. We will use units h = 1, except for 
the current defined in Eq. ([5]). 

The two conduction bands can be classified into the 
s-wave and p-wave parts, defined by 



,J R C e,R* + v L C eXa 



J L^e,Ra 



J R^eXa 



R^ V l 



(4) 



The s-wave electrons couple to the dot via the tunneling 
matrix elements, as Ht — J2a ^ { d t^ua- + ^ta d cr)' wnere 
S„ = J de^J p{e) and v 



v L . The p-wave 



electrons do not couple to the dot, but they contribute 
to the current through the impurity^ 



J 



_ v R J L +v 2 L J R 



4^E(^-4^)- (5) 



Here, P a = J de\J p(e) p ea , and Jl (Jr.) is the current 
flowing from the left lead (dot) to the dot (right lead). 
The scattering states in the s-wave part, cv eo ., can be 
calculated explicitly for Hjj — including all effects of 
Ht, 



= S' 



v^pJe)G r (e)4 
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id 



(6) 
(7) 



Note that the resonance width of the impurity level is 
given by A(e) = nv 2 p(e) =T L + T R with Y v = irvlp(e) 
in the noninteracting case. We will also be using the 
advanced Green's function Gq(s) = {G (e)}* . The one- 
particle states a ea and p £a . constitute the eigenstates of 

Hq = H Q + Hq , 

Hi = ^2jdee a\ a a ea , Hg = ^ J, deep\ a p e<y . (8) 

cr cr 

The scattering state is also normalized as {a e a , a\, ,} = 
$ era 1 5(e — s'), and the impurity state d a can be expressed 
in terms of a Ea as 



d„ =v dey/p(e)G r (e) a E 



(9) 



Here, we have omitted, for simplicity, the possibility that 
bound states could emerge outside of the conduction 
band. It can be taken into account when it is necessary. 



B. Nonequilibrium steady state 

The noninteracting Hamiltonian can also be diagonal- 
ized with the scattering states moving towards left and 
right, 



H °= dse lt,uale,u<r i 

T n ~ J 



v=L,R cr 



le,Ra 



V R a Ety + v L p ea 



2 T n ,2 ' 'e 



%,L* 



(10) 

V L a ea- v RPea 



(11) 



This propagating-wave form of Hq are compatible with 
the standing- wave form Eq. JHJ), as the scattering states 
have two-fold degeneracy. The chemical potentials pl 
and p R can be assigned to these two propagating states 
lt,L<r and it, Ren respectively, using 1 ^^ 



^0 = ^" j d£ ilt,Lo-le,Lo ~ lt,Ra%,Ra] 

+ tJt^| E/ & ( a - a -~^-)- (12) 

With this operator, the noninteracting density matrix 
can be constructed as 



Po 



E a = Tie-^ H "- Y ^ . (13) 



The statistical average is defined by (• • • )o = Tr \p ■ ■ ■ ]. 
For instance, the distribution function for 7 6 va is given 
by 

(lLo-%'v'o->)o = U{e) S(e - e') 8 VV , 8 aa> , (14) 
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fes(e) 



FIG. 1: Distribution function / c ff(e) at T = 0. 



where /^(e) = f(e — p u ), and /(e) = [e £ / T + l] is the 
Fermi function. Similarly, from Eqs. pip and (|14[) . we 
obtain 



("Low)o =/eff(e)5(£-e')^a 



/eff(e) 



_ r £ / £ ( £ ) + r fl / fl (e) 
r L + r R 



(15) 
(16) 



and 



[WeJ-^l^-fi')^- (17) 



Here, the distribution function f e s{s) appears not only 
for the impurity state but also for all the scattering states 
with the s-wave character £ At zero temperature, f c s{s) 
takes a constant value Tl/^r + Tl) for p R < e < fiL, 
as shown in Fig. [TJ Specifically, the constant becomes 
1/2 for the symmetric coupling Tl = Tr. This feature 
looks quite similar to a high-temperature behavior of the 
Fermi function /(s)^^ = 1/2, and it can be explained 
as a manifestation of the mixed-state nature of the RDM, 
which is given in Eq. ([33]). 

The Coulomb interaction is switched on adiabatically 
based on the Keldysh formalism^ 



p = 5(0, -co) p S(- oo,0) , 
r t 2 

S(t 2 ,ti) = Tcxp 



dte^Hue-^* 



—I I aid "nut 
ti 



(18) 
(19) 



Note that [5(0, — oo), p ea \ = 0, and thus the interaction 
affects only the s-wave part of the Hilbert space, which 
is described by the Hamiltonian 



(20) 



Hershfield has provided one possible way to construct an 
interacting operator Y, by which the full density matrix 
is written in the form p oc e ~P( H - Y ) , 23 In the present 
study, however, we go back to the original definition of p 
given in Eq. (|18[) for a general description. 

The interacting average (• • • ) is taken with the full den- 
sity matrix p. For instance, using Eq. ([5]), the average 



for the local charge can be written in the form 
(n d ) = v 2 J2 [dsdeWp(e')p(e)G Q a (e')G° r (e){at, a a Ea ), 



Here, G<(ui) = i f^dt e iut (dl e' 



H b t, 



iH b t 



) is the 



lesser Green's function, which for the noninteracting 
limit is given by G<(e) = -f cS (e) [G r {e) - Gg(e)]. The 
nonequilibrium current can also be expressed as an aver- 
age defined in the s-wave subspace^ 



(J) 



G r (e) (p| v a £CT ) - Gg(e) (a^p^) 
* Jdu[f L (u)-f R (u)]^^-[-TmCr{u;)]. 



(22) 



Here, G r (u) = -i j^dte^+^^e^d^-^^dt}) 
is the full retarded Green's function. Therefore, to cal- 
culate these averages, the p-wave degrees of freedom can 
be integrated out first before taking the trace over the 
correlated s-wave part. 



III. REDUCED DENSITY MATRIX 

We describe in this section a derivation of the RDM. 
To this end, we introduce a discretized Hamiltonian cor- 
responding to Eq. {SJ as Hq = Hq + Hq, 



a m 

^0 =Zm £ ™ (pt m ,aPe m ,<T - (P. 



9 S \ (23) 



,o-/0 J 



(24) 



r>} = im'L', and 



The constant term (• • • )[f s , 



with the normalization {ct Sm CT , a £ 

{Pe m ,a>Pl m , ,cr>} = <W'<W 

which denotes the noninteracting ground-state average 
in equilibrium, is subtracted so as to make the lowest 
energy of Hq and that of Hq zero, 



\ GS - (rf n ) GS 

,o70 — \Pe m ,<rPe m ,a/0 



1, £m < 0, 

0, e m > 0. 



(25) 



Specifically, we assume a linear discretization with a fi- 
nite level spacing Se, such that e m — (m + 1/2) Se for 
772 — 0, ±1, ±2, . . .. The continuous spectrum can be re- 
covered in the limit of Se — > with a E a = 5 £m a / VSe 

and p e a = p e a j\fSe. The discretized operator is also 
introduced for the propagating states 7 e va using Eq. 
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(fTTjl . and then Hq can also be written in the form 

Ho = EEE e ™ (jL t utr% m ,ua - <7L,^%,„^a>o :s 



v a m 



(26) 



In this discretization, each discrete level e m preserves the 
two-fold degeneracy due to the left-going and right-going 
scattering states. This degeneracy is essential to describe 
the current-carrying state with the density matrix. 

The trace over the Hilbert space can be carried out 
explicitly for the discretized model. The normalization 

factor So = Tr e^ /3( - H °~ Y °\ where Yq is the discretized 
version of Eq. (fT2|l , can be calculated as 



— o — 



n nn[ 

u=L,R iy e m 



(27) 



Here, we have used a property e^^ L+ ^ R ^ = 1 following 
from = — fiji = eV/2. The corresponding free energy 
that is defined by Sq = e~@ F ° takes the form 



f3F 



1 

"Se 



E E h ln 



u=L,R a 



-f}(e—H„) sgne 



(28) 



In the continuum limit Se — * 0, the factor l/5e should be 
replaced by the Dirac ^-function 5(0). It appears because 
the energy of the whole system is proportional to the 
system size N. For instance, the level spacing is given 
by Se = 2D/N for a flat band with a half width D. 
Therefore the level spacing Se has the physical meanings 
for the averages of order N, although we have started 
from the continuous model in order to capture properly 
the two-fold degeneracy in the scattering states. 
The density matrix can also be calculated as, 



,-/3(H -Y ) 



Po 



T± e -P(H -Y ) 

nnn 

v=L,R a e 



1 + e-Pie-Vv) 



(29) 
(30) 



Here, we have suppressed the suffix m of the discrete fre- 
quency e m , for simplicity. Equation (|30[) can be rewritten 
in terms of a e a and p s a , by substituting Eq. (JTTJ), as 



?o=nn( i -^)( i -/«) 

a e 



l 



-Ik !l_ + /?_ _ i ) 5 t 5 



h 



h 



vl + v\l-f L vl + vll-f R 



l)pl 



PeaVea 



VLVR 



!l - fR 



vl + v* R (l-f L )(l-f R ) 
(l-2/ L )(l-2/ fl ) ^ _ ^ _ 



(31) 




FIG. 2: (color online). Parameter /3 e fi = 1/T e , which deter- 
mines the RDM as Eq. (|34[) , is plotted as a function of e for 
several temperatures T/(eV) = 0.0, 0.05, 0.1 and 0.2, in the 
symmetric coupling case Fl = Tr. 



Furthermore, the partial trace over the p-wave part can 
be carried out to yield the RDM for the s-wave subspace 



Trpp =nn( i -/° ff ) 



1 



2.f, 



eff 



1 



q ct 

Lt e<r L scr 



1 - /off 

It can also be expressed in an exponential form 



(32) 



Po = Trp po = 



Trc e K o 



K " = EE t b ^ E / de | 

(7 £ (T 

l-/eff(e)" 



■ln 



/eff(e) 



(33) 

:<r> 

(34) 
(35) 



Then, the average of an operator Os = Os(a, a'), 
which is an arbitrary function of a £ a and a\ a , can be 
calculated in the s-wave subspace using the RDM, 



(0 s )o = Tr s Tr p p O s 



= Tr c 



(36) 



Note that the impurity level has been included in the 
s-wave part. It is straight forward to confirm that 
{a\ a a e , CT ,)o given in Eq. (p~5|) can be reproduced from 
Eq. (|3"6")) . and then the function / e ff(e) defined in Eq. 
(|16p can be expressed in terms of T E as 



/off(e) = 



1 



+ 1 



(37) 



The bias dependence appears in the s-wave subspace 
through this energy-dependent parameter T E . The in- 
verse of this parameter, /3 e ff = l/T e , is plotted in Fig. [2] 
as a function of e for several values of real temperature 
T in the symmetric coupling = Tr case. The energy 
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dependence of T e becomes weak with increasing T . At 
zero temperature, it takes the form, 



% 



oo, |e| < eV/2, 
0, |e| > eV/2. 



(38) 



Namely, T e reaches infinity inside the energy window 
|e| < eV/2, while it remain zero on the outside. For 
this reason, in the limit of eV — > oo the nonequilib- 
rium Green's functions at the impurity site coincide 
with the equilibrium ones for the high-temperature limit 
T — > oo.— 

The partial trace over the p-wave part can be carried 
out also for Hjj ^ 0, as the Coulomb interaction affects 
only the s-wave part of the Hilbert space. Therefore, 
the interacting RDM p s evolves adiabatically from the 
noninteracting RDM as 



= S , (0,-oo)^5(-oo,0) , 



(0 s )=Tr s Tr p 



pO s 



Tr c 



(39) 
(40) 



It enables us to treat the bias voltage within the s-wave 
part of the Hilbert space. Namely, one can start with 
the single-channel Anderson model H s in the continu- 
ous form, Eq. (|20[) . Then, the noncquilibrium distribu- 
tion can be taken into account through the initial RDM 
jojf , which evolves to the interacting one p s . For the 
single-channel model, it is not necessary to use the same 
discretization scheme that was introduced to obtain Eq. 
(|32| . because pjf has a well-defined continuum form as 
Eq. (f34|) . In addition to these points, the correlation func- 
tions containing a p-wave operator can also be calculated 
with the single-channel model as shown in Appendix [A"l 



IV. MIXED-STATE ASPECT FOR T L = T R 

As we see in Fig. [2] effects of the bias voltage are 
pronounced significantly at zero temperature, especially 
when the system has an inversion symmetry Tl = Tr. 
We consider in the following the properties of the RDM 
in this particular case. 

At zero temperature T = the noninteracting density 
matrix p for the whole system can be constructed as a 
pure state I'&ncq), m which the scattering states incident 
from left and right are filled up to p,^ = eV/2 and p,R = 
—eV/2, respectively, 



i<*w=n n iHl^^o) , (4i) 
i*o> = n ek^v io> - n n syt i°> • ^ 



e<0 a 



£<0 cr 



|0) is the vacuum state. The noninteracting ground 
|$ ) can be decomposed into the Fermi sea of the 
s-wave part |$q) an d that of the p-wave part l^o") as 



Here, 
state 



|$ ) = | Therefore, Eq. (gU can be rewritten, 

using Eq. (TiTj) for v L = vr, 



as 



— 

n n n 



(SL -Pea) ( a e'a+Pe'a) 



V2 



V2 



iE'l«f>l«0 • 



l*o) 
(43) 



Here, |$f ) and |$f ) are the excited states of the s-wave 
part and that of the p-wave part, respectively: the sign 
caused by exchanges of the fermion operators is absorbed 
into the global phase of each many-particle eigenstate in 
the definition. The primed sum for i runs over a set 
W, which consists of all the possible excited states that 
can be generated from |$ ) by changing the occupation 
of the one-particle states inside the biased energy region 
|e| < eV/2. The factor Zq corresponds to the dimension 
of W. 

The pure-state average of an arbitrary operator Os = 
0s (a, o^) that is defined in the s-wave space can be 
calculated, carrying out the summation over the p-wave 
part, as 



($ncq|Os 



ijeW 



= ^Y,'( (p ?\ d s\^)= TT s\p S oVs}- (44) 



tew 



Therefore in this case, the noninteracting RDM takes a 
mixed-state form Pq = Y^ n \^n) w o,n(^n\ w ^ n 



u 0,n 



o, 



neW, 



(45) 



The summation with respect to n, or Tr s [- • •], run over 
the whole region of the s-wave subspace. The statistical 
weight WQ n takes just two values or l/Z . However, 
it is not a simple function of the many-particle energy 
eigenvalue E^ n of H$ |*£ > = |^>. 

There is another important relation between $f) and 
| $f ) , as these states are generated simultaneously in Eq. 
(|43f . If a one-particle level e in the biased energy region 
|e| < eV/2 is being occupied by an s-wave electron a\ a 
in |$f ), then the corresponding level has to be unoccu- 
pied in the p-wave part \<&f ), or vice versa. Therefore, 
there exists a pair of states, i and i' , that show a relation 
««,|*f>|*f) = -p^l*!,)^) for | e | < eV/2. Here, 
the minus sign is caused by the order of the operators 
for the one-particle level e. Using this property, a p-wave 
operator appearing in an average can be replaced by an 
s-wave operator as 

(<J>aeq| C S ?W|$„c q ) = ^ £ '<*f |($f | 6 S p £(7 \^pf) 



ijew 



= - lf L (e) - f R (e)} Tr s p s Q O s 



(46) 



6 



Similar relation holds also for the asymmetric coupling 
Tl 7^ Tfl and at finite temperatures as shown in Ap- 
pendix El 

The interacting density matrix at T — is also con- 
structed as a pure state with the wavefunction 




(47) 

Here, = 5(0, — oo) ), and it is an interacting 

eigenstate^ (H$ + H v ) = E% ). The wave 
function for the p-wave part keeps the noninteracting 
form, because 5(0, -co) commutes with the p-wave op- 
erators. Therefore, the summation over the p-wave part 
can be carried out in the same way as that in Eq. (|4"4")l . 
and the interacting average for the operator Os can be 
written in the form, 

(* neq ia s |* neq ) = ^^'(*f|a s |*f) 

Z/n * — ' 

= J2 w L«\6s\0- (48) 

n 

The interacting RDM in this case also takes a simple 
form p s = J2 n \^n) w o,n(^n\- Equation (JHJ gives us 
much insight into the steady state. It will provide, ob- 
viously, an interacting equilibrium average if we replace 

w o,n by Peq n = e ~^ En I^Z n e ~ l3En ■ Therefore, some of 
the nonequilibrium properties could be inferred from the 
difference between the statistical weights Wq u and yof qjn . 

As Eq. (|4"5)) is being expressed in terms of the interact- 
ing eigenstates, Hamiltonian-based nonpcrturbative ap- 
proaches developed for equilibrium, such as exact diag- 
onalization and NRG, can be applied to the matrix el- 
ement (Nf^lOg \^n)' Furthermore, the bias voltage ap- 
pears just through Wg n , which keeps the noninteracting 
form given in Eq. (|45p . Thus, what is necessary to carry 
out the summation with respect to n in Eq. (|4"5)) is the 
information about whether or not each interacting eigen- 
state belongs to the set W in the limit of Hu — > 0. 
Namely, the correspondence between |^f) and its non- 
interacting counterpart \$f) is an additional issue that 
one has to solve in the nonequilibrium case. At low en- 
ergies, the correspondence can be treated correctly with 
the quasi-particles of a local Fermi liquidi 7 - ' 13 ! 14 ' 16 ! 17 ' 18 

The evolution of the wavefunctions is described for- 
mally with the time-dependent operator 5(0, —oo). Al- 
ternatively, it can be expressed in a time-independent 
form as that in a stationary-state scattering theory ; 26 i 27 
and several versions of the time-independent approaches 
have been examined with some different points of 
viewAi— — — — In our description, the adiabatic evo- 
lution of the interacting eigenstates from the noninter- 
acting ones could also be traced directly in the Hilbert 
space of a discretized version of the single-channel model, 
by changing the coupling constant U continuously. 
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FIG. 3: (color online). Averaged statistical weight Wq(E) = 
W(E) I Q(E) vs many-particle excitation energy E, where 
N = eV/5e. 



V. ENERGY DISTRIBUTION IN THE MIXED 
STATES 

We consider in the following the mixed-sate properties 
for the symmetric coupling = Tr further at T = 0. 
The pure state |3> ne q) given in Eq. (|TTj) is an eigenstate, 
Ho |$ncq) = E Q lcq I'&ncq), of the noninteracting Hamilto- 
nian with an excitation energy 

eV/2 ( v\ 2 

E Mi^rw- < 4<J) 

e m =-eV/2 

Here, we retained only the dominant term that is pro- 
portional to the system size of order TV. As mentioned in 
Sec. QUI the level spacing is given by Se = 2D /N for a flat 
band with a half width D. The feature of the pure state 
reflects in the s-wave subspace through the RDM. The 
mean excitation energy in the subspace is equal to a half 

of the one in Eq. ([35]). Eq = Tr s [p£#jf] = E^ cq /2. Note 
that the contribution of the impurity level is of order 1, 
namely a 1/N portion of the total energy. 

In order to clarify how the nonequilibrium weight w$ n 
distributes in the many-particle Hilbert space consisting 
of the eigenstates we examine an averaged weight 

Wq(E) over an equienergy surface; 

Wq (E) = W(E)/Cl(E) , (50) 

W(E) =J2' S ( E ~ O. n (E) = £ S ( E - O (51) 

ieW n 

It is normalized as Zq = JdE wg(E)Sl(E), and may be 
used as (* ncq |a s |* ncq ) «£„^M(*n|Os|*2)/Z 
for a rough approximation. We calculate wJq (E) for sev- 
eral finite 6e, and the results are shown in Fig. [3] the 
parameter Nq = eV/Se corresponds to the number of 
one-particle states in the energy window |ej < eV/2, and 
thus Z — 4 N ° . The averaged weight Wq (E) becomes 
a constant for small E, as all the many-particle excited 
states for < E < eV/2 belong to the set W defined 
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in Eq. (|43[) . This flat ^-dependence is quite different 
from the exponential dependence of the Boltzmann fac- 
tor e~ /3E in the thermal equilibrium. We also see that 
Wq (E) penetrates into the higher energy region over the 

bias energy eV. It spreads out, approximately, up to 
s 

E ~ cE Q with a numerical constant c of order 10 _1 . 

Note that the mean excitation energy is proportional to 
s 

N as E — N eV/4. We have also obtained the analytic 
expressions of Sl(E) and W(E) for small Se as shown in 
Appendix [B] and have checked that the same behavior 
could be seen for larger iVo assuming a Gaussian form for 
W(E) given in Eq. (|B7[) . instead of the exact integration 
form Eq. |B6|) . 

The feature of Wq (E) seen in the s-wave subspace is 
caused by the fact that the bias voltage has been in- 
troduced not into the many-particle energy but into the 
one-particle states as the width of the energy window 
between the two chemical potentials pl and /xr. It is 
quite different from the way the temperature determines 
the equilibrium distribution e~ E / T , in which T scales di- 
rectly the many-particle energy E. The original weight 
Wq n must have the similar features in the n dependence. 
However, whether or not the high-energy many-particle 
states give significant contributions to the nonequilib- 
rium averages depends also on the properties of the op- 
erator Os- For instance, the noninteracting correlation 
functions, such as (af a a e , CT ,)o and Gq(e), can be re- 
produced from Eq. (|4"4")) . The feature of w^ n reflects in 
these correlation functions just through / c fr(e). It should 
be regarded, however, those are the results that are de- 
termined by all the many-particle states belonging to the 
set W. Therefore, it is necessary for calculating the inter- 
acting Green's functions with Eq. (j4"5)) to take all these 
high-energy states into account. The Kondo resonance 
at finite bias voltages is being affected by these features 
of the noncquilibrium distribution. 



We have studied further the properties of the RDM at 
zero temperature for the symmetric coupling Tl = Tr. 
In this case the noncquilibrium averages can be written 
in a Lchmann-representation form Eq. (|48p with a sim- 
ple mixed-state distribution wfn given in Eq. (|45[) . It 
has also been deduced that the mixed states contain large 
portions of high-energy many-particle states. Our formu- 
lation could be used for constructing Hamiltonian-based 
nonperturbative approaches in future. 
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APPENDIX A: RDM APPROACH TO THE 
p-WAVE OPERATOR 

The correlation functions which include a p-wave op- 
erator can be expressed, using Eq. (|3ip , as 



{£>SPea)o = lim Tr c Tr, 



v l v r fL(s)-f R (e) 



Tr, 



., 5 FnOsa ea (Al) 

J R T v L JeS{ £ ) 1 J 

in the noninteracting case. For instance, Eq. (|17[) can be 
reproduced from Eq. (|A1[) for Os = o<^ E , al . The interact- 
ing average can also be expressed in a similar form 



(O sPsa ) = -- 



v l v r 



/eff(e) 



x Tr s pg S(-oo, 0) Os 5(0, -oo) a £ 



(A2) 



VI. SUMMARY 

In conclusion, we have reexamined the nonequilibrium 
steady state for a quantum dot, checking out the sta- 
tistical distribution in the many-particle Hilbert space. 
We have, first of all, discretized the conduction bands 
keeping the two- fold degeneracy due to the left-going and 
right-going scattering states. This degeneracy is essen- 
tial to describe the current-carrying state driven by the 
bias voltage. Then, by integrating out the p-wave part of 
the channel degrees of freedom, the noninteracting RDM 
Pq has been obtained explicitly as shown in Eq. f|33|) . 
It brings the information about the bias voltages to the 
remaining s-wave subspace, and determines the nonequi- 
librium distribution. The interacting RDM p s , which is 
defined in Eq. (|39p . evolves adiabatically from p^. Using 
the RDM, the steady-state averages Eq. (l4"0"|) can be cal- 
culated, in principle, with the single-channel Anderson 
model. 



APPENDIX B: n(E) AND W(E) FOR SMALL Se 

We provide the integration forms of £l(E) and W(E) 
for small level spacing Se. To this end, we consider a 
grand partition function Q{t, E c ) for noninteracting elec- 
trons in a flat band with the cut-off E c . The grand 
partition function and many-particle density of states 
F(E, E c ) are relating each other with the Laplace trans- 
form, 

j* oo 

G(t,E c )= dE e- TE F(E,E C ) . (Bl) 
Jo 

For small Se, the left-hand side, G(t, E c ), is given by 

a A=e,/i£=0 a 

(B2) 

X(r,E c )^[ E d eln(l + e-n = m ~ HTEc) . (B3) 
Jo T 
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Here. 



. (*) = E™ =1 ^-e— and 0(0) = vr 2 /12. 
The label A represents the contributions of the electron 
(e) and hole (h) excitations. The function F(E, E c ) is 
given as the inverse Laplace transform 



-tE c 



Tq—IOO 



2vri 



) 2 4 0(O)-0(rg c ) 
>— e Ji r , 

(B4) 



The function {1(E) = F(E, oo) is obtained taking the 

0, 4>(tE c ) -> 0, and the 



limit of E r 



then 



-tE c 



integration can be carried out for E > 0, 



1 40(0) 



E5e 



(2^jA(j){Q)E/5i 



(Wo) 



Here, I\(x) is the modified Bessel function, which for 
large x takes the asymptotic form; Ix(x) ~ e x /^/2ttx. 

The other function W(E) = F(E, eV/2) is obtained 
for E c = eV/2. Note that W(E) = {1(E) for < E < E c 



because the cut-off does not affects the excitations in 
this energy region. Furthermore, W(E) = for E > 
2E C (N C + 1) beyond an upper bound, where N c = E c /Se. 
Equation (|B4[) can be rewritten taking the path along 
the imaginary axis as 



iEt 



1 



-UE C 



4 4>{Q)-<f>{itE c ) 

-e3i ~t . (B6) 



Note that Re [0(0) - <j>(ixj\ = x 2 /4 for — 7T < X < 7T. 
Thus, W(E) has a peak at E = E C (N C + 1), and around 
the peak it may be approximated in a Gaussian form 



W(E) 



*2N a 



E r 



n(2N c + 3) 



e ( 2N - 



^t(^-^i) 2 



(B7) 



This approximate form overestimates W(E) for small E, 
suggesting that the real density given by Eq. (|B6|) has 
larger weight around the peak than that of the Gaussian. 
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